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ABSTRACT 

Ongoing and planned wide-area surveys at optical and infrared wavelengths should detect a few times 10 s galaxy 
clusters, roughly 10% of which are expected to be at redshifts > 0.8. We investigate what can be learned about 
the X-ray emission of these clusters from the ROSAT All-Sky Survey. While individual clusters at redshifts > 0.5 
contribute at most a few photons to the survey, a significant measurement of the mean flux of cluster subsamples 
can be obtained by stacking cluster fields. We show that the mean X-ray luminosity of clusters with mass M>2x 
1O 14 /z~'M0 selected from the Sloan Digital Sky Survey should be measurable out to redshift unity with signal-to- 
noise > 10, even if clusters are binned with Az = 0. 1 and AlnM ~ 0.3. For such bins, suitably chosen hardness ratio 
allows the mean temperature of clusters to be determined out to z ~ 0.7 with a relative accuracy of AT /T < 0.15 
forM> 10 14 /t -1 M ro . 



1. Introduction 
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With moderately deep, wide-area imaging surveys in the opti- 
cal or near infrared it is now possible to detect large samples of 
galaxy clusters. Dalcanton (1996) proposed that clusters could 
be detected as surface brightness enhancements even when all 
but a few of their galaxies are too faint to be detected indi- 
^ vidually. Her suggested procedure consists of identifying and 
*^ removing stars and galaxies from carefully flat-fielded images, 
smoothing the residual with a kernel matched to the core size 
O-of clusters, and searching for significant peaks in the resulting 
q smoothed map. Gonzalez et al. (2001) successfully constructed 
5_i the Las Campanas Distant Cluster Survey (LCDCS) by applying 
this technique to drift-scan data taken with the Las Campanas 
(7$ Great Circle Camera (Zaritsky, Schectman & Bredthauer 1996). 
; * They mapped well over 100 square degrees and constructed a 
. catalog of 1073 groups and clusters. The estimated redshift lim- 
k> its of the catalog range from ~ 0.3 for groups to ~ 0.9 for mas- 
j_j sive galaxy clusters. 

d ' The high intrinsic uniformity of drift-scan surveys like the 
LCDCS makes them ideal for applying Dalcanton's cluster- 
detection technique. In a theoretical study, Bartelmann & White 
(2002) showed that massive galaxy clusters should be detectable 
in the Sloan Digital Sky Survey (SDSS) out to redshifts of ~ 1 .2 
if data in the r 1 , i' and z' bands are summed. For the final pro- 
jected SDSS survey area of 10 4 square degrees, > 10 5 galaxy 
clusters should be detectable at the 5-G level, and ~ 10% of 
those are expected to be at redshifts > 0.8. 

Relatively little is known about the X-ray emission of clus- 
ters at redshifts beyond ~ 0.5 despite numerous cluster surveys 
based on X-ray data. The main reason for this is the steep de- 
crease with redshift of the observed X-ray flux, which implies 
that at z > 0.5 individual massive clusters produce at most a 
few photons in surveys like the ROSAT All-Sky Survey (RASS; 
Snowden & Schmitt 1990). Detections of distant galaxy clusters 
from the X-ray data alone are limited to very luminous systems 
which can be detected in the restricted areas where deeper ob- 
servations are available. 

The upcoming availability of large cluster surveys in wave- 
bands other than the X-ray regime allows a reversal of the tradi- 
tional X-ray survey strategy. Rather than identifying clusters in 
the X-ray data, it becomes possible to stack X-ray survey data 
for a large number of fields where clusters are already known 
from other surveys. The low background count rate at X-ray 



wavelengths makes this an efficient technique for detecting the 
summed emission from a large stack of clusters. 

In this paper we investigate the prospects for using the RASS 
to detect X-rays from suitable samples of clusters identified in 
the SDSS data. In Sect. ^]we describe our model for the cluster 
population. This is based closely on the properties of nearby 
clusters and specifies the cluster distribution in mass, redshift, 
X-ray temperature and luminosity. In Sect. ^] we convert cluster 
X-ray luminosities to count distributions expected in the ROSAT 
All-Sky Survey. Based on these, we calculate in Sect. K\ the 
expected signal-to-noise both for the detection of mean cluster 
emission and for estimates of mean cluster temperature. Sect. |5] 
summarises and discusses our conclusions. 



2. Model specifications 

2.1. Cosmology 

Much evidence suggests that the universe is spatially flat with 
low nonrelativistic matter density £2o- Baryons make up only a 
small fraction of this matter; the rest is dark, presumably con- 
sisting of some massive, weakly interacting particle. A cosmo- 
logical constant £2a or an equivalent "quintessence" field con- 
tributes the remaining energy density. For definiteness, we as- 
sume £2o = 0.3, £2a = 0.7 and h = 0.7 

We assume structure to form from an initially Gaussian den- 
sity fluctuation field 8 with statistical properties specified by 
its linear power spectrum, for which we adopt the CDM form 
given by Bardeen et al. (1986) with primordial spectral index 
n = 1. The only remaining free parameter is then the normali- 
sation of the initial fluctuation field which we take as Os = 0.9. 
This value was originally estimated based on the observed lo- 
cal abundance of galaxy clusters (White, Efstathiou & Frenk 
1993; Eke, Cole & Frenk 1996; Viana & Liddle 1996; Pierpaoli, 
Scott & White 2001; Evrard et al 2002) but some recent anal- 
yses favour smaller values (Reiprich & Bohringer 2002; Viana, 
Nichol & Liddle 2002; Lahav et al. 2002). We will show results 
for (jg= 0.9 ±0.1. 



2.2. Cluster population 

Haloes form from Gaussian primordial density fluctuations 
through gravitational collapse. Press & Schechter (1974) first 
derived an approximate formula for the mass distribution of 



haloes as a function of redshift z. This has recently been modi- 
fied by Sheth, Mo & Tormen (2001) and Sheth & Tormen (2002) 
based on an ellipsoidal rather than a spherical model for col- 
lapse. They give the differential comoving number density of 
haloes as 



2 ( 1 \ 6 dv ( v l 
n(M,z)6M =A\ - 1 + ^5- — — exp ) AM . 



(1) 



where v = ^/abcO^ 1 (M,z) defines the linear amplitude required 
for collapse of a density fluctuation and p is the mean cosmic 
density today. a(M,z) in this definition is equal to Gq(M)D + (z), 
where Oo(M) is the present rms fluctuation in the dark matter 
density contrast within spheres containing the mean mass M, 
and D + (z) (with D + (Q) — 1) is the growth factor for the lin- 
ear growing mode (cf. Carroll, Press & Turner 1992). The lin- 
ear density contrast required for collapse 8 C depends weakly on 
cosmology; for the ACDM model we have chosen 8 C = 1.673 
(e.g. Lokas & Hoffman 2001). The parameters A, a and q are 
constants; the original Press-Schechter formula is obtained from 
([!]) by putting A = 0.5, a = 1 and q = 0. This mass function, 
with A = 0.322, a — 0.707 and q = 0.3, has been shown to fit 
high resolution numerical simulations of structure growth in a 
wide range of cosmologies, provided the halo mass is defined at 
fixed density contrast relative to the cosmic mean density (Jenk- 
ins et al. 2001). 

Next, we need to know the X-ray luminosity of a cluster of 
mass M. We adopt the observed relation between cluster tem- 
perature T and bolometric X-ray luminosity Lx 



Lx = 10 44 ergs- 1 



kT \ 
1.66keV J 



2.331 



(2) 



as derived by Allen & Fabian (1998). Observations suggest that 
there is little evolution in the Lx — T relation out to redshifts 
z ~ 0.4 (e.g. Mushotzky & Scharf 1997; Allen & Fabian 1998; 
Reichart, Castander & Nichol 1999). Lacking any reliable in- 
formation about evolution to higher redshifts, we assume @) to 
hold at all redshifts. This, of course, is a major uncertainty of 
our study. 

According to the virial theorem, a halo of mass M in equilib- 
rium at redshift z with a structure similar to observed clusters 
should have a mean temperature given by 



ytr=4.88keV 



Mh(z) 



W 5 Mr. 



2/3 



(3) 



where h(z) is the Hubble constant at redshift z in units of 
lOOkms -1 Mpc -1 and, in contrast to Eq. ([!]), M is here defined 
as the mass interior to a sphere with mean overdensity 200 times 
the critical value at redshift z. Recall that we assume h(0) = 0.7 
throughout our analysis. The constant in this relation is taken 
from the cluster simulations of Mathiesen & Evrard (2001 ; their 
Table 1) and is appropriate for specifying the temperature of the 
best fit single temperature model for the X-ray spectrum over the 
mass and redshift ranges of interest. When necessary, we use an 
NFW model of concentration parameter 5 to convert between 
cluster masses defined at different overdensities. 



2.3. X-ray emission 

We assume that clusters emit X-rays through thermal radia- 
tion according to a Raymond-Smith plasma model (Raymond 
& Smith 1977). Apart from cluster temperature and redshift, 
the model has two free parameters, the metal abundance and 



an overall normalisation corresponding to the total X-ray lu- 
minosity. We fix the metal abundance to Z = 0.3 Z Q at all z 
in agreement with the observed abundances of local clusters 
(e.g. Fukazawa et al. 1998). The results of Schindler (1999) sug- 
gest little evolution towards higher redshift and the final count 
rates we derive depend only very weakly on metallicity. Thus 
neglecting any dependence on redshift does not induce signifi- 
cant uncertainty. 

Let F v (T,z)dv be the total X-ray luminosity emitted in the 
spectral interval [v, V + dv] by a cluster of temperature T at red- 
shift z. If the cluster is observed in an energy band bounded 
by Ei and E% >E\, only a fraction / of its bolometric flux is 
included in the bandpass, where 



f= / F v (T,z)dv 



F v (T,z)dv 



-l 



(4) 



Thus the band-limited flux Sx is related to the bolometric X-ray 
luminosity through 



Sx = 



fLx 
4%Dl(z) 



(5) 



where £>l(z) is the luminosity distance from the observer to red- 
shift z. Note that this flux must still be modified to account for 
foreground absorption. 

We use version 11.1 of the xspec software package (Arnaud 
1996) to tabulate / for an observing band between 0.5 and 
2.4keV, for cluster temperatures between 0.5 and 12keV, and 
for redshifts between and 2. Interpolating within this table and 
using Eqs. (|2|), (|j) and (|1), we can convert cluster masses to clus- 
ter temperatures, X-ray luminosities, and finally to unabsorbed 
fluxes in the observed energy range. 

The azimuthally averaged X-ray surface brightness profile 
E(9) of galaxy clusters is often modelled using the so-called beta 
profile (Cavaliere & Fusco-Femiano 1978), 



E(9) = E 



-(3P-1/2) 



(6) 



where 9 C is an angular core radius, and the amplitude Eo is cho- 
sen to produce the required X-ray flux Sx- Based on observation, 
we choose fi = 2/3 (e.g. Mohr, Mathiesen & Evrard 1999). For 
the linear core radius r c , we adopt the relation 



125kpc/; 



-l 



Lx 



5 x 10 44 ergs- 



0.2 



(7) 



where Lx is the X-ray luminosity between 0.5 and 2.4keV. This 
relation is a fair representation of at least some clusters with lu- 
minosities within 10 43 ~ 45 ergs~' (Jones et al. 1998). Following 
Vikhlinin et al. (1998), we assume that (Q) does not evolve with 
redshift. The angular core radius is then 9 C = r c D~ l (z), where 
D(z) is the angular-diameter distance. In fact, Eq. (^) is a poor 
fit to the profiles of many clusters, particularly those with strong 
apparent cooling flows. This is not, however, of any great conse- 
quence for our modelling since the RASS does not, in any case, 
resolve the inner regions of most clusters. 

Having fixed p\ Sx and the angular core radius 9 C , the beta 
profile is normalised by 



E = 



Sx 

27T92 



(8) 



3. Halo detection 

3.1. Point-spread function 

The point-spread function /(6;£,<|>) of the ROSAT-PSPC had 
three components, a Gaussian kernel, Lorentzian wings, and a 
component which falls off exponentially with the angular sepa- 
ration 8 from the centre of the image (Hasinger et al. 1995). The 
parameters for these components generally depend not only on 
photon energy E, but also on (j), the off-axis angle of the source. 

The width of the PSPC point-spread function can be charac- 
terised by the effective solid angle S£2 (E,§) covered, 

8Q.(E,ty)=2n[ 9d9/(9,£,<M , (9) 
Jo 

and we can define an effective radius 9 e ff (£,0) by 

M£ ,, ) = («)"\ (10) 

The effective radii for six different off-axis angles between 10' 
and 60' are shown as functions of photon energy in Fig. [l]. 
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FIG. 2. — The ROSAT-PSPC point-spread function, averaged over off- 
axis angles within the PSPC field-of-view. The different curves show 
the point-spread function for four different photon energies, as indi- 
cated. 



an average effective radius valid for the hard band of the All- 
Sky Survey. Performing this average and weighting the photon 
energies with the effective detector area as a function of E, we 
find 9 eff = 2.1'. 
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FIG. 1 . — Effective radii of the PSPC point-spread function as functions 
of photon energy, for off-axis angles between 10' and 60' (from bottom 
to top). 

The field-of-view of the PSPC was large, with a radius of ap- 
proximately 60'. Since a given point on the sky was scanned 
at many different off-axis angles during the All-Sky Survey, the 
appropriate point-spread function for the ROSAT All-Sky Sur- 
vey at a given photon energy is an area-weighted average of 
f(B;E,§) over the field-of-view, 



ME) = 



(60') 2 Jo 



60' 



4>d4>/(e,£,4>) 



(ii) 



Figure || shows the result for four different photon energies be- 
tween 0.5 and 2.0keV. 

Figure || shows that the point-spread function, averaged over 
off-axis angles, falls to ~ 10% of its peak value within 2-3 arc 
minutes with little dependence on photon energy. Figure [l] con- 
firms this weak dependence on photon energy, and illustrates 
the strong dependence of effective PSF radius on off-axis angle. 
While 9 e ff is below 0.5' for nearly on-axis photons, it increases 
above 1' for photons coming from the edge of the field-of-view. 

The effective radius of the averaged point-spread function 
f(Q,E) can finally be averaged over photon energies to obtain 



3.2. Converting fluxes to count rates 

We now need to estimate the signal expected in the All-Sky Sur- 
vey from a cluster with unabsorbed flux Sx and temperature T at 
redshift z. To do this, we first modify the fluxes Sx calculated us- 
ing Eq. (||) to allow for absorption by foreground neutral hydro- 
gen. We assume a constant hydrogen column of 4 x 10 20 cm -2 , 
which is typical for the high galactic latitudes covered by the 
SDSS (e.g. Dickey & Lockman 1990). We convert the absorbed 
fluxes to PSPC count rates, using the fakeit task of the xspec 
package with the PSPC response matrix]]. 

In practice, we compute a two-dimensional table containing 
PSPC count rates in the energy range between 0.5 and 2keV 
for absorbed Raymond-Smith spectra of a fixed unabsorbed 
flux normalisation and for cluster temperatures between 0.5 and 
12keV and redshifts between and 2. Fluxes determined from 
(|j) can then be converted to absorbed count rates by interpolat- 
ing within this table. 

3.3. Exposure times; background level 

The effective exposure time in the All-Sky Survey varies across 
the sky because of the ROSAT scanning strategy. It is high- 
est near the ecliptic poles and lowest close to the ecliptic plane 
(cf. Snowden et al. 1995). Maps for the exposure time and the 
background count rates were downloaded from the ROSAT All- 
Sky Survey web page[| The left panel in Fig. || shows the cumu- 
lative exposure-time distribution for the complete All-Sky Sur- 
vey (dashed curve), and for the area around the Northern Galac- 
tic cap covered by the Sloan Digital Sky Survey. The median 
exposure times are marked by vertical lines. 

The effective exposure times on the whole sphere and on the 
SDSS area are only marginally different. For the SDSS area, we 
find a median value 

fexp=414s. (12) 



electronically provided at ftp://ftp.xray.mpe.mpg.de/rosat/calibration/data 
2 http://www.xray.mpe. mpg.de/cgi-bin/rosat/rosat-survey/ 




FIG. 3. — Cumulative distributions of exposure-time (left panel) and 
background count-rate (right panel) in the ROS AT All-Sky Survey. The 
dashed curves refer to the complete survey, the solid curves to the SDSS 
area only. The vertical lines mark the medians. 



Similarly, the background level is anisotropic across the sky. 
The right panel in Fig. || shows the cumulative distributions of 
the background count rate in the All-Sky Survey for the whole 
sky (dashed line) and for the SDSS area (solid line). 

The background count rate within the SDSS area is noticeably 
lower than on the whole sky; its median value is 



B = 0.94s" 1 deg- 



2.61 x lO-VVcmin- 2 



Results 



(13) 



Figure |] shows photon-count contours in the plane spanned by 
cluster mass and redshift. On a grid covering that plane, we com- 
pute temperature, luminosity, flux, and count rate as described 
in the previous section. We then multiply the count rate by the 
median exposure time in the SDSS area, averaged over photon 
energies in the 0.5 — 2.4keV band. 

The contours are logarithmically spaced by 0.25 dex between 
0. 1 and 100 counts (upper and lower solid curves, respectively). 
They appear jagged because a substantial fraction of the X-ray 
flux is contributed by metal lines which move in and out of the 
observed energy band as the redshift changes. The contours be- 
come smooth if the metal abundance is set to zero. From this 
plot one can see, for example, that the redshift limit below which 
individual clusters contribute more than ten photons to the All- 
Sky Survey increases from z m ax ~ 0.1 at M = 10 14 hr x M Q to 
Zmax ~ 0.8 at M — 1O 15 /i _1 M . The dashed curve shows the 
redshift limit for detection of clusters as 5-G surface brightness 
enhancements in the combined r 1 , i 1 and z' data of the SDSS 
(Bartelmann & White 2002). Clearly the SDSS should produce 
cluster catalogues which are much deeper at all masses than 
those that can be made from the RASS. 

Figure Q illustrates that only 0.3 photons per cluster are ex- 
pected for clusters ofM~ 1O 14 A _1 M at redshift z ~ 0.8. The 
number of such clusters expected in the SDSS is so large, how- 
ever, that it should be possible to determine their mean X-ray 
properties by stacking data for many fields. This is true even 
if the mass-redshift plane is divided into relatively narrow bins. 
We now investigate this in more detail. 

The background level of the All-Sky Survey is quite low, of 
order ls~'deg~ 2 which translates to approximately 0.8 total 
counts per resolution element within the median exposure time 



FIG. 4. — Contours in the mass-redshift plane showing the counts re- 
ceived per cluster within the effective radius of the PSPC point-spread 
function. The contours are logarithmically spaced at 0.25 dex between 
0.1 (upper solid contour) and 100 counts (lower solid contour). The 
dashed curve marks the expected upper redshift limit for 5-CT cluster 
detection in the combined r 1 , i' and z' bands of the SDSS. The con- 
tours for low-mass clusters appear jagged because their X-ray spectra 
have strong features due to heavy elements which move relative to the 
observed energy band. 



of the survey. The background will nevertheless dominate the 
noise in a stacked image of distant clusters. Let B be the mean 
surface density of background photons in a single image, and 
C(M,z) be the expected number of photons from a single cluster 
of mass M at redshift z. Let p(Q) be the expected surface den- 
sity of these cluster photons as a function of angular distance 
from cluster centre. p(Q) is given by a convolution of the 
mean cluster surface brightness profile [Eq. (||)] with the point- 
spread-function of the survey (Fig. |2|) and we normalise it so 
that / p(0)27t9d8 = 1. In practice for distant clusters the p.s.f. 
is much broader than the cluster image so that p(Q) is propor- 
tional to the p.s.f. itself. 

For a stack of N cluster fields the surface density of the back- 
ground is NB and the expected surface density profile is NCp (8) . 
Assuming Poisson photon statistics, the optimal estimator of the 
cluster signal is then: 



NC= / w(Q)[0(d)-B]2nQdd , 



(14) 



where 27i(3(9)9d9 is the observed photon count in an annu- 
lus width d9, and the filter function w, normalised so that 
/ w(d)p{B)2%edB = 1, is given by 



w(8) 



p(e) 



p(d)+B/C 



p 2 2nQdQ 



p + B/C 



-l 



(15) 



Clearly the expectation value of the estimator of equation ([U 
is just NC while its variance is 



Var(AfC) =NC / w 2 (9)/?(9)27r9d9 . 



(16) 



Thus the expected signal-to-noise for detecting the stacked clus- 
ter is 



/7 2 (9)27i9d9 



P (B)+B/C 



1/2 



(17) 



If clusters are individually well above background (Cp(Q) ^> 
B over most of the broadened image) this gives the obvious 
result, (S/N) « (NC) l l 2 for the stack. When background 
dominates (Cp(O) <C B) the corresponding result is (S/N) « 
NC/{NB J p 2 (e)27i9de] 1 / 2 . In both cases the signal-to-noise of 
the detection grows as N l l 2 for the stacked image. Figure [| 
shows the number of cluster fields required for a 5-a detection 
in the stacked image as a function of cluster mass and redshift. 



function (hi) multiplied by the comoving cosmic volume, 
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FIG. 5. — Number of clusters to be stacked to achieve a 5-0 detection of 
their total X-ray emission in the ROSAT All-Sky Survey. The contour 
levels are 1 (lower contour), 10, and 100 (upper contour). The dashed 
curve marks the expected upper redshift limit for 5-G cluster detection 
in the combined r', i' and z' bands of the SDSS. 

Contours are shown for N = 1 (lower solid curve), TV = 10, 
and N = 100 (upper solid curve). The figure shows that it 
takes 100 stacked cluster fields to obtain a 5-a detection of 
clusters with M — 10 14 ft -1 M Q at redshift z ~ 0.4, but the con- 
tours rise steeply enough that with the same number of stacked 
fields one reaches redshifts above unity for cluster masses M > 
3 x 1O 14 /i _1 M . As in Fig. |], the dashed line shows the upper 
redshift limit expected for 5-a cluster detection in the combined 
r 1 , i' and z' bands of the SDSS. 

We now have to compare the number of cluster fields needed 
to achieve a high signal-to-noise ratio with All-Sky Survey data 
to the number of clusters we can expect to be available. The idea 
is to select fields around clusters which are known from other 
data, and we continue to take the SDSS as an example. We 
therefore ask how many clusters can be expected in the SDSS 
data. 

To give specific examples, we select two redshift intervals of 
width Az — 0.1 each, one over 0.6 < z < 0.7 and the other over 
0.9 < z < 1.0. Our previous work has obtained the expected 
redshift limit zn m (M) as a function of cluster mass for detec- 
tion in SDSS data (Bartelmann & White 2002). For a 5-C de- 
tection in the combined r J , i' and z' bands, it is indicated by 
a dashed line in Figs. [I] and |j. For each redshift interval, we 
thus know the completeness limit in cluster mass, i.e. the low- 
est cluster mass M\\ m above which clusters in that interval are 
expected to be detectable. For the lower and upper redshift in- 
tervals defined above, we obtain Mr nn = 3.9 x 10 13 /! _1 M Q and 
Mn m = 2.0 x 10 A -1 M0 respectively. For each interval we then 
define a series of mass bins between Mr lm and 10 15 h~ l M & such 
that AlnM~ 0.3. 

The number of clusters in the redshift interval [zi,Zi + Az] per 
mass bin [Mj,Mj + \] is obtained through an integral of the mass 



H+Az rM j+ i 



ANij= / dz / dMn(M,z) 



dV 



dz 



The volume per unit redshift is 



dV 



dz 



■ %D 2 (z) 



dD 



prop 



dz 



(1+z) 3 . (18) 



(19) 



where D is the angular diameter distance and D prop the proper 
distance. The factor K instead of 4k accounts for the fact that the 
SDSS only covers a quarter of the sky. Figure || shows the re- 
sulting cluster numbers AN,j and the total photon numbers AC; 
expected from these clusters. The solid and dotted curves show 
results for the lower and upper redshift bins, respectively. In 
order to illustrate the sensitivity of the results to Os, the error 
bars mark the range obtained for Gg — 0.9 ±0.1. The curves 
showing the total photon counts received in each mass bin are 
flatter than those showing the total cluster number because clus- 
ters with higher mass are more X-ray luminous. 
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FIG. 6. — Number of clusters AN (upper panel) and total cluster counts 
AC (lower panel) in the two redshift intervals 0.6 < z < 0.7 (solid curve) 
and 0.9 < z < 1.0 (dotted curve) in mass bins of logarithmic width 
AlnM ~ 0.3 between the SDSS completeness limit in the respective 
redshift interval and 10 i5 A _1 Mq. The total counts received from all 
clusters per mass bin drop much less steeply than the cluster number 
because the number of counts received per cluster increases strongly 
with cluster mass. The error bars bracket results obtained by changing 
Og by ±0. 1 and illustrate the sensitivity to the power- spectrum normal- 
isation. 

The figure shows that, even with relatively fine mass binning, 
more than 10 4 clusters should be detectable per mass bin be- 
low 1O 14 /! _1 M in the lower redshift interval 0.6 < z < 0.7. 
For comparison, Fig. |] shows that several hundred stacked clus- 
ter fields are already sufficient for a 5-a X-ray detection in the 
RASS. Similarly, more than 10 3 clusters are expected per mass 
bin below 4 x 1O 14 A -1 M at higher redshifts, 0.9 < z < 1.0, 
where fewer than ~ 100 cluster fields need to be stacked for an 
X-ray detection. A useful way to quantify these numbers is by 
calculating the expected signal-to-noise ratio in a stack of all the 
cluster fields in each mass bin and in each of our two redshift 
intervals. The results are shown in Fig. ^| 

At the lower redshift the signal-to-noise ratio starts above 40 
near 4 x 10 13 /i _1 Mp,, where the contribution of metal lines to 
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FIG. 7. — Signal-to-noise ratios in stacked cluster fields in the given 
mass bins for the two redshift intervals 0.6 < z < 0.7 (solid curve) and 
0.9 < z < 1.0 (dotted curve). As in Fig. ^, the error bars show the 
range obtained by varying a$ by ±0.1. The signal-to-noise ratio in the 
lower redshift interval reaches ~ 40 near 3 x l0 u h~ 1 M Q . Even near 
10 15 ft _1 Af© in the upper redshift interval, the signal-to-noise ratio is 
~ 10. 



the flux is high. With increasing mass, the line contribution de- 
creases and S/N has a shallow minimum near 10 14 A -1 Af Q . In- 
creasing continuum emission causes a broad peak at > 40 cen- 
tred on 3 x 10 14 /z _1 M Q . It then decreases slowly towards higher 
masses. The drop-off results from from the low cluster number 
at the high-mass end. If we set the metal abundance to zero, the 
low X-ray flux at the low-mass end makes the signal-to-noise 
drop to ~ 20 near 4 x 10 13 Ii~ 1 Mq. Even in the upper redshift 
interval, the signal-to-noise ratio is above 10, rising to > 20 in 
the lowest mass bin. These results are, however, very sensitive to 
G8. Near lO 15 /z _1 M in the upper redshift interval, the signal- 
to-noise ratio varies between ~ 5 and ~ 20 as Og is increased 
from 0.8 to 1.0. 

The high signal-to-noise ratio even for high-redshift clusters 
encourages us to investigate whether it will be possible to esti- 
mate cluster temperatures from hardness ratios. We introduce 
two energy bands, one with 0.5 < E /keV < 1 and the second 
with 1 < E /keV < 2. The counts C\2 in these two bands deter- 
mine the hardness ratio 



hard counts 
soft counts 



Ci 
Ci 



(20) 



We use xspec to compute the hardness ratio l^(T,z) expected 
for RASS data for clusters with temperature T at redshift z- For 
clusters of mass M at redshift z, the uncertainty in the tempera- 
ture measurement is then 



AT{M,z) 



9g. 
dT 



T(M),z] A* 



-l 



(21) 



where the uncertainty A^, of the measured hardness ratio (|20j) is 
determined by the count statistics. The boundaries of the energy 
bands were chosen so that 3i is typically of order unity in the 
mass and redshift ranges considered here. The signal-to-noise 
ratio of the hardness ratio ^,/A^, is > 10 for all cluster mass 
bins in the redshift interval 0.6 < z < 0.7, and is > 8 for the bins 
in the redshift interval 0.9 < z < 1.0. The derivative of %. with 
respect to T is ~ 0.8 for T ~ IkeV and falls to ~ 0.1 the high- 
est temperatures. As a result temperature determinations should 



be most accurate for clusters with M ~ lO 14 /z _1 M ; at lower 
masses, line emission in the low-energy band dominates and the 
uncertainty A$t increases because of poor photon statistics in the 
high-energy band. We show T and AT/T in Fig. || for the same 
mass bins and redshift intervals used previously. For compari- 
son, the plot also gives the mean cluster temperature expected 
as a function of mass in each redshift interval. Note that both T 
and AT/T are emitted rather than observed values. 



0.5 - 



0.1 




0.05 



FIG. 8. — The curves with open squares show the relative uncertainty 
AT/T of cluster temperatures determined from the hardness ratio be- 
tween a soft (E e [0.5, l]keV) andahard(£ e [l,2]keV) band. Clusters 
are stacked in mass bins in the two redshift intervals 0.6 < z < 0.7 (solid 
curve) and 0.9 < z < 1.0 (dotted curve). As in Figs, pi and ffl, the error 
bars indicate the range obtained by varying o"g by 3:0.1. \Vhile the 
temperature uncertainty is very large for the low-mass clusters, it drops 
near 10% for moderate-redshift clusters withM > 1O 14 A _1 M0, and is 
< 20% for the high-redshift clusters with masses < 6 x 10 14 h~ l M G . 
The curves with filled circles show the cluster temperature in keV for 
the given mass bins and redshift intervals. 

Figure [| shows that the relative uncertainty in the mean tem- 
perature of the clusters in each mass bin is remarkably small for 
0.6 < z < 0.7. For cluster masses > 10 14 /i _1 M it is AT/T < 
0.15. Over the mass range 10 14 — lO 15 /z _1 M , it appears that 
a > 10a measurement of cluster temperature should be possi- 
ble. As in Figs. || and 0, the error bars in Fig. |^ indicate the 
range obtained by varying the power-spectrum normalisation c?8 
by ±0.1. For clusters in the high-redshift band, 0.9 < z < 1.0, 
the relative temperature uncertainty increases both because of 
count statistics and because of decreasing sensitivity of to T . 
Despite this, temperature measurements at 3 to lOo should be 
possible. Note that a careful maximum likelihood measurement 
of T would give results with somewhat higher significance than 
the simple hardness ratio approach we have adopted here. 



5. Discussion 

Ongoing and planned wide-area surveys will detect tens of thou- 
sands of galaxy clusters out to redshifts near and above unity. 
For example, searching for surface-brightness enhancements in 
a smoothed stack of the r 1 -, i'- and z'-band data of the Sloan 
Digital Sky Survey should allow clusters of 5 x lO 13 /z _1 M to 
be detected out to z ~ 0.7, while z > 1 is reached for masses 
above - 3 x lO 14 /z _1 M (Bartelmann & White 2002). 

We have investigated here whether existing X-ray data can be 
used to measure the X-ray emission of these clusters by stack- 
ing sufficiently many fields. We assume clusters to be distributed 



in mass and redshift as given by the numerical results of Jenk- 
ins et al. (2001). Their temperatures are taken to be propor- 
tional to [M h(z)] 2 ^ , with the normalisation taken from the N- 
body/SPH simulations of Mathiesen & Evrard (2001). We adopt 
the observed low -redshift relation between bolometric X-ray lu- 
minosity and temperature, and we assume that it holds at all 
redshifts. We model cluster X-ray surface-brightness profiles by 
a beta profile, although this has little effect on our results be- 
cause most distant clusters are not resolved in the RASS. We 
convert the bolometric X-ray luminosity into a count rate using 
the xspec software, assuming a Raymond-Smith plasma model 
with a metallicity of 0.3 solar and a foreground neutral-hydrogen 
column of 4 x 10 20 cm~ 2 . 

The only suitable survey of the X-ray sky is the ROSAT All- 
Sky Survey (RASS). With its median exposure time of approxi- 
mately 415 seconds and its effective detector area of ~ 230 cm 2 , 
it detected ~ 10 photons from a cluster of mass 10 14 A -1 Af s at 
z ~ 0.1, and only about one photon from a similar cluster at 
z ~ 0.5. Since the effective angular resolution of the RASS is 
~ 2', cluster emission is typically spread over an effective solid 
angle of ~ 14 square arcminutes. Due to the low background 
of the PSPC detector, only ~ 1 .5 background photons are ex- 
pected within this solid angle during the median RASS exposure 
time. This corresponds to the number of photons expected from 
a cluster with mass M~ 10 14 A -1 Af a at redshift z~ 0.35, or with 
mass M~4x 10 u h~~ l M @ at z ~ 1. Thus stacked cluster fields 
are background dominated at lower mass or higher redshift than 
this. 

Requiring a signal-to-noise exceeding 5, we find that 100 
fields must be stacked to get a significant detection of clusters 
withM~ \Q u hr l M @ at z ~ 0.4, or with M ~ 3 x 10 14 A -1 Af Q 
at z ~ 1 -0. A stack of ten cluster fields should give a 5-G detec- 
tion of massive clusters with M ~ 1 1 5 h~~ l M @ at z ~ 1.1, should 
any such exist. 

The number of clusters expected in wide-field surveys like the 
SDSS is enormous and allows the detection of X-ray emission 
from even fairly low-mass clusters out to surprisingly high red- 
shift. In the redshift interval between 0.6 and 0.7, the surface- 
brightness technique of Dalcanton (1996) should detect clusters 
in the SDSS data down to a mass limit of — 3.9 x 1O 13 A -1 M . 
If we bin the clusters by mass into logarithmic bins with width 
A In M = 0.3, the signal-to-noise ratio for the X-ray detection ex- 
ceeds 35 near rises above 40 near 3 x 10 14 /z _1 M Q 
and drops to ~ 25 at 10 /i -1 M©. In the interval between red- 
shifts 0.9 and 1.0, the mass completeness limit for SDSS cluster 
detection increases to 2 x 10 ft M©, but X-ray detections are 
still possible with signal-to-noise ratios above 10 if clusters are 
binned by mass as described. 

We have also shown that the signal-to-noise ratio of the 
stacked cluster images is high enough to divide the photons into 
two energy bands, E E [0.5,1] and E e [1,2], and to estimate 
cluster temperatures from the count ratio. In particular, for clus- 
ters with 0.6 < z < 0.7 and masses > 10 14 /! _1 M Q , the hardness 
ratio changes with cluster temperature sufficiently strongly for 
mean cluster temperatures to be determined with a typical rela- 
tive uncertainty of AT /T < 15%. 

Of course, these results depend on the modelling assump- 
tions we have made. Their sensitivity to changes in the power- 
spectrum normalisation Og is shown in Figs. ||, ^ and 0, where 
the error bars bracket results obtained adopting Og = 0.9 ±0.1. 
Other critical assumptions are that the relation between bolomet- 
ric X-ray luminosity and temperature is independent of redshift, 
and that the cluster temperature scales with cluster mass as given 
by simulations. Our assumptions about cluster X-ray profiles are 
less critical because of the low angular resolution of the RASS. 



Using photometric redshifts for brightest cluster members, it 
should be possible to determine redshifts for SDSS clusters with 
an accuracy of Az ~ 0.05. On the other hand, at best very rough 
estimates of cluster mass can be obtained from the optical data. 
The optically selected clusters in a given redshift interval can 
be binned by magnitude, and the study suggested here will then 
give relations between optical luminosity and mean X-ray lumi- 
nosity and temperature. The latter can then be used to give an 
improved estimate of mean cluster mass. 

Wide-area surveys in the microwave regime will be carried 
out in the near future which will detect of order one cluster 
per square degree trough the thermal Sunyaev-Zel'dovich effect. 
The Planck satellite, for instance, due for launch in early 2007, is 
expected to detect of order 30000 galaxy clusters on the sky out- 
side the Galactic plane, approximately 10% of which will be at 
redshifts beyond 0.5. Stacking these clusters in the same way as 
described here, and combining their total integrated Compton-y 
parameter with their X-ray emission, will allow their total bary- 
onic mass and perhaps their temperatures to be constrained. 
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